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Theoretical approaches that use one-body densities as dynamical vari- 
ables, such as Hartree-Fock or the density functional theory (DFT), break 
isospin symmetry both explicitly, by virtue of charge-dependent interac- 
tions, and spontaneously. To restore the spontaneously broken isospin 
symmetry, we implemented the isospin-projection scheme on top of the 
Skyrme-DFT approach. This development allows for consistent treatment 
of isospin mixing in both ground and exited nuclear states. In this study, 
we apply this method to evaluate the isospin impurities in ground states of 
even-even and odd-odd N ~ Z nuclei. By including simultaneous isospin 
and angular-momentum projection, we compute the isospin-breaking cor- 
rections to the + — > + superallowed /3-decay. 
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1. Introduction 

The atomic nucleus is a quantum system composed of the two types of 
strongly interacting fermions, the nucleons. The charge independence of 
the nuclear interaction is at the roots of the isospin symmetry [Tj [2] . This 
concept remains valid even in the presence of Coulomb interaction, which is 
the major source of the isospin breaking. This is so because of the smallness 
of the isospin-breaking isovector and isotensor components of the Coulomb 
field as compared to the isospin-conserving components of the nuclear and 
Coulomb forces. 

The isotopic spin quantum number, T, provides strong selection rules 
for nuclear reactions, decays, and transitions [3]. In particular, the selection 
rules for /3-decay Fermi and Gamow- Teller transitions are AT = and AT = 
0, ±1, respectively, with the exception ofT = 0— > T = transitions that 
are forbidden [H [5]. The superallowed + — > + Fermi transitions bridge 
nuclear structure with the electroweak standard model of particle physics, 
providing the most accurate estimate for the V u d matrix element of the CKM 
matrix hence, testing the CKM unitarity. From a nuclear structure 

perspective, the unitarity test depends critically on a set of theoretically 
calculated isospin-breaking corrections whose precise determination poses a 
challenging problem [HJ El [TO] . 

In this work, we calculate the isospin impurities and isospin-breaking 
corrections to the superallowed Fermi decay by using a newly developed 
isospin- and angular- momentum-projected DFT approach without pairing 
[11[ [T2| I13J . This technique takes advantage of the ability of the mean 
field (MF) to properly describe long-range polarization effects. The MF 
treatment is followed by the isospin projection to remove the unwanted 
spontaneous isospin mixing within MF [12\ W$\ [T5"l I16|. 

This paper is organized as follows. We begin in Sec. [2] with a short sum- 
mary of our isospin- and angular- momentum-projected DFT approach. In 
Sec. El we present applications of the isospin-projected DFT variant of the 
model to the isospin mixing in the ground states (g.s.) of even-even N=Z 
nuclei. Section H] discusses preliminary results for the isospin-breaking cor- 
rections to the superallowed beta decays calculated by considering simulta- 
neous isospin and angular-momentum restoration. Finally, the conclusions 
are contained in Sec. 

2. Theory 

The isospin-projected DFT technique [HI [121 E3] utilizes the ability of 
the self-consistent MF method to properly describe the balance between 
the long-range Coulomb force and the short-range nuclear interaction, rep- 
resented in this work by the Skyrme-type energy density functional (EDF). 
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To remove the spurious isospin-symmetry-breaking effects, we use the stan- 
dard one-dimensional isospin projection after variation, which allows us to 
decompose the Slater determinant |$) into good isospin states \T,T Z ): 

m= E b T,T z \T,T z ), £ |6t,tJ 2 = 1. (1) 

T>\T Z \ T>|T Z | 

Here, Py r stands for the conventional one-dimensional isospin-projection 
operator: 

\TT Z ) = -J—Rf m 



2^/Ntt z Jo 



r d(3 T sin/3 T dr, T ,0*r) A(#r)l*)> (2) 



where /3t denotes the Euler angle associated with the rotation operator 
R{Pt) = e~ l ^ TTy about the y-axis in the isospace, d?p T {Pt) is the Wigner 
function [17J, and T z = {N — Z)/2 is the third component of the total isospin 
T. The normalization factors Ntt z , or interchangeably the expansion coef- 
ficients &t,t z that encode the isospin content of |<£), read: 

N TTz = |6r,rJ 2 = ($|P^ T J$) 

2T + i r T 



r d(3 T sin fa 4 zTz (fir) A/"(/3 T ), (3) 
Jo 



10 

where AT(Pt) = {*&\R(fir)\$) is the so-called overlap kernel. For technical 
aspects concerning the calculation of the overlap and Hamiltonian kernels, 
we refer the reader to Ref. [13]. The isospin-projected DFT technique uti- 
lizes the ability of the HF solver HFODD [18] to produce fully symmetry- 
unrestricted Slater determinants [$). 

The isospin projection determines the set of good isospin states (called 
the basis in the following) , which in the next step is used to rediagonalize the 
entire nuclear Hamiltonian, consisting of the kinetic energy, Skyrme EDF, 
and the isospin-breaking Coulomb force. The rediagonalization leads to the 
eigenstates: 

\n,T z )= ^t,tJT,T z ), (4) 

T>\T Z \ 

numbered by index n. The amplitudes aj- r z define the degree of isospin mix- 
ing through the so-called isospin-mixing coefficients (or isospin impurities) 
for the n— th eigenstate: 

«c = 1 " K.tJLx, (5) 



4 
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where l^yy (max stands for the squared norm of the dominant amplitude in 
the wave function \n,T z ). It is worth stressing that the isospin projection, 
unlike particle-number or angular-momentum projections, is essentially non- 
singular; hence, it can be safely used with the local EDFs. The rigorous 
analytical proof of this useful property can be found in Ref. [13] . 

The combined isospin and angular-momentum projection leads to the 
set of states, 

\I,M,K;T,T Z ) = 1 PlrPLxm, (6) 

\JNtT z ;IMK 

which form another normalized basis built on |3>). Here, Pj> t ano - ^MK 
stand for the isospin and angular-momentum projection operators, respec- 
tively and M and K denote the angular-momentum components along the 
laboratory and intrinsic z-axes, respectively |19j . Now the problem becomes 
more complicated because of the overcompleteness of the basis §6§ related 
to the X-mixing. This is overcome by performing the rediagonalization of 
the Hamiltonian in the so-called collective space, spanned for each / and T 
by the natural states, \IM;TT z )( l \ as described in Refs. [181 EO]- Such a 
rediagonalization gives the solutions: 

\n;IM;T z )= ]T a™\IM; 7T Z >«, (7) 
i,T>\T z \ 

which are labeled by the index n and by the conserved quantum numbers 
I, M, and T z = (N - Z)/2 [cf. Eq. ©]. 



3. Isospin mixing 

By using the perturbation theory [22] and the analytically solvable hy- 
drodynamical model [23] , the isospin mixing in atomic nuclei has been stud- 
ied since the 1960s (see Ref. [24J for a review). These simple approaches 
accounted for such qualitative features of the isospin impurities like the 
steady increase in N = Z nuclei with increasing proton number and strong 
quenching with increasing \N — Z\. Quantitatively however, their predic- 
tions for the values of the isospin impurities ac were not very reliable. 

Increased demand for accurate values of isospin mixing has been stim- 
ulated by the recent high-precision measurements of superallowed /3-decay 
rates [61 [7] . Large-scale shell- model approaches [25J , although very accu- 
rate in the description of configuration mixing, can hardly account for the 
long-range polarization exerted on the neutron and proton states by the 
Coulomb force whose accurate treatment requires using large configuration 
spaces. In contrast, in self-consistent DFT, such polarization effects are 
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Fig. 1. Isospin impurities in the 
ground state of 80 Zr, predicted 
by DFT, using various Skyrme 
parametrizations [H| plotted versus 
the corresponding excitation ener- 
gies of the T = 1 doorway states. 
Open dots mark results obtained be- 
fore the Coulomb rediagonalization 



(BR), a^ H) = 1 - |& r=|r ,| J Tj a > 



29.5 30.0 30.5 31.0 31.5 



Ex=! [MeV] 



which were calculated by using ex- 
pansion coefficients of Eq. ([T]). Full 
dots mark the impurities ([5]) ob- 
tained after the Coulomb rediago- 
nalization (AR). 



naturally accounted for by finding the proper balance between the Coulomb 
force, which tends to make the proton and neutron states different, and the 
isoscalar part of the strong force, which has an opposite tendency. 

In general, isospin impurities determined without removing spurious 
isospin mixing are underestimated by about 30% compared to the values 
obtained after rediagonalization [12J. In the particular case of 80 Zr, the 
removal of spurious admixtures increases ac from ~2.9% to ~4.4%, as il- 
lustrated in Fig. [TJ It is encouraging to see that the latter value agrees 
well with the central value of empirical impurity deduced from the giant 
dipole resonance 7-decay studies, as communicated during this meeting by 
F. Camera et al. [26J. Unfortunately, experimental error bars are too large 
to discriminate between various Skyrme parametrizations, which differ in 
predicted values of ac by as much as ~10%. 

Figure [2] illustrates our attempts to correlate the values of ac with the 
surface and volume symmetry energies, which are primary quantities char- 
acterizing the isovector parts of nuclear EDFs. The linear regression coeffi- 
cients shown in the figure hardly indicate any correlation of ac with these 
quantities. In fact, no clear correlation was found between the calculated 
values of ac and other bulk characteristics of the Skyrme EDFs, including 
the isovector and isoscalar effective masses, and incompressibility. 

4. Isospin-breaking corrections to the Fermi matrix elements of 

the superallowed /3-decay 

An accurate evaluation of ac is a prerequisite for determining the iso- 
spin-breaking correction 5c to the + — > + Fermi matrix element of the 



S U 0.8 
0.7 



SkO 



SIII 



^ o 

SLy 

R= 0.29 



SkO 



sni 



R= 0.68 



18 19 20 21 22 
a sym (/W2) [MeV] 



18 19 20 21 22 
a sym (/W2) [MeV] 



1.1 

^ 10 
^ 0.9 

S U 0.8 

0.7 



40 Ca 



SkP 



MSkl 
^.SkM* 



SLy 



SLy5 



StOQr 



SkO'% 



sm 



SkO 

R= 0.08 



100 Sn 



SkP# 



MSkl 



SIII 



Skkc 



SkM* 



SLy SLy5 



SkO" 



sfo 
R=0.10 



29 30 31 32 
a sym (/W [ MeV l 



29 30 31 32 
a sym (/W [MeV] 



Fig. 2. Isospin impurities predicted by several Skyrme EDFs for the ground states of 
40 Ca (left) and 100 Sn (right) plotted versus the surface (top) and volume (bottom) 
symmetry energy. 



isospin raising/lowering operator T± between nuclear states connected by 
the superallowed /3-decay: 

\{F = 0+,T « 1, T z = ±l\f ± \r = 0+, T « 1, T z = 0)| 2 = 2(1 - S c ). (8) 

Here, the state \I W = + ,T & 1,T Z = ±1) corresponds to the g.s. of the 
even-even nucleus whereas \I W = + ,T ~ 1,T X = 0) denotes its isospin- 
analogue in the neighboring N = Z odd-odd nucleus. Unlike the former 
one, the odd-odd configuration cannot be expressed in a form of a MF 
product wave function [13]. Therefore, to compute the states in odd-odd 
N = Z nuclei, we use the following strategy (see Fig. 1 of Ref. [13J for a 
schematic illustration) : 

• Firstly, we compute the so-called antialigned g.s. configuration, |i/(g)7r) 
or |f®7f), by placing the odd neutron and the odd proton in the lowest 
available time-reversed (or signature-reversed) single-particle Nilsson 
orbits. 

• Secondly, to correct for the fact that the antialigned configurations 
manifestly break the isospin symmetry, that is, |z^(g)7r) ~ — = 0} + 
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Fig. 3. Isospin impurities in Sc, cal- 
culated for four antialigned configura- 
tions that are obtained by putting the 
valence neutron and proton in opposite- 
K Nilsson orbitals originating from the 
/V/2 shell, that is, \vK ® irK) with 
K = 1/2,3/2,5/2, and 7/2. Open and 
full dots show the results obtained by 
employing only the isospin projection 
and simultaneous isospin and angular- 
momentum (/ = 0) projection, respec- 
tively. 
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Fig. 4. Isospin-breaking correction to 
the Fermi matrix element for the super- 
allowed transition 14 0— J> 14 N. Full dots 
represent our results plotted as a func- 
tion of the basis size (number of HO 
shells taken in HF calculations) . A con- 
servative 10% error was assigned to the 
last two points. The values quoted in 
Refs. [7] (including errors) and [27] are 
shown for comparison. 



\T = 1)), we apply the simultaneous isospin and angular- momentum 
projection to create the good isospin and good angular momentum 
basis \I, M, K, T, T z = 0) of Eq. ©. 

• Finally, to obtain the state |J = 0, T ~ 1,T Z = 0), we rediagonalize 
the total Hamiltonian, including the Coulomb term, in the new basis 
[cf. Eq. 47}]. 

The projected |P r = + ,T ~ 1,T Z = ±1) states in even-even nuclei are 
computed in the same way. 

Restoration of angular momentum turns out to be the key ingredient 
in evaluation of the isospin impurity in odd-odd nuclei. This is illustrated 
in Fig. [3l which shows ac calculated for the T ~ 1 states in 42 Sc. Four 
solutions shown in Fig. [3] correspond to the four possible antialigned MF 
configurations built on the Nilsson orbits originating from the spherical 
vfl/2 and vr/ 7 / 2 subshells. These configurations can be labeled in terms 
of the K quantum numbers, K = 1/2, 3/2, 5/2, and 7/2, as \vK <g> irK). 
In a simple shell-model picture, each of those MF states contains all /=0, 
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2, 4, and 6 components. Prom the results shown in Fig. [33 it is evident 
that the isospin projection alone (upper panel) leads to unphysically large 
impurities, whereas the impurities obtained after the isospin and angular- 
momentum (1=0) projection (lower panel with the scale expanded by the 
factor of 500) are essentially independent of the initial MF configuration, 
as expected. The average value and standard deviation of 0.586(2)% shown 
in the figure were obtained for the configuration space of N = 10 spherical 
harmonic oscillator (HO) shells, whereas for N = 12 the analogous result is 
0.620(2)% (see below). 

Although indispensable, the angular-momentum projection creates nu- 
merous practical difficulties when applied in the context of DFT, that is, 
with energy functional rather than Hamiltonian. The major problem is 
the presence of singularities in energy kernels [25]. Although appropriate 
regularization schemes have already been proposed [25], they have neither 
been tested nor implemented. This fact narrows the applicability of the 
model only to those EDF parametrizations which strictly correspond to 
an interaction, wherefore the singularities do not appear. For Skyrme-type 
functionals, this leaves only one EDF parametrization, namely SV [30] • This 
specific EDF contains no density dependence and, after including all ten- 
sor terms in both time-even and time-odd channels, it can be related to a 
two-body interaction. Despite the fact that for basic observables and char- 
acteristics such as binding energies, level densities, and symmetry energy, 
SV performs poorly, we have decided to use it in our systematic calculations 
of 5c- Indeed, while SV would not be our first choice for nuclear structure 
predictions, it is still expected to capture essential polarization effects due to 
the self-consistent balance between the long-range Coulomb and short-range 
nuclear forces. 

In order to test the performance of our model, we have selected the su- 
perallowed /3-decay transition 14 — )- 14 N. This case is particularly simple, 
because (i) the participating nuclei are spherical and almost doubly magic, 
which implies suppressed pairing correlations, and (ii) the antialigned con- 
figuration in 14 N involves a single \vp\/2 ® npi/2) configuration that is 
uniquely defined. The predicted values of 5c are shown in Fig. [J] as a func- 
tion of the assumed configuration space (that is, the number of spherical 
HO shells N used). While the full convergence has not yet been achieved, 
this result, taken together with other tests performed for heavier nuclei, 
suggests that at least N = 10 shells are needed for light nuclei (A < 40), 
whereas at least N = 12 shells are required for heavier nuclei. The resulting 
systematic error due to basis cut-off is estimated at the level of ~10%. 

Even though calculations for all heavy (A > 40) nuclei of interest are 
yet to be completed, and due to the shape-coexistence effects there are still 
some ambiguities concerning the choice of global minima, our very prelim- 
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inary results are encouraging. Namely, the mean value of the structure- 
independent statistical-rate function Ft, obtained for 12 out of 13 tran- 
sitions known empirically with high precision (excluding 38 K— 7- 38 Ar case), 
equals Ft = 3069.4(10), which gives the V ud = 0.97463(24) amplitude of the 
CKM matrix. These values match very well those obtained by Towner and 
Hardy in their latest compilation [7]. That said, owing to the poor quality 
of the SV parameterization, the confidence level [10J of our results is low. 
On a positive note, our method is quantum mechanically consistent (see 
discussion in Ref. [5]) and contains no free parameters. 

5. Summary 

In summary, the isospin- and angular-momentum-projected DFT the- 
ory has been employed to calculate isospin mixing and isospin-breaking 
corrections to the + — > + Fermi superallowed /3-decay. Our parameter- 
free model capitalizes on the ability of the MF approach to describe long- 
range polarization effects. The self-consistent HF wave functions containing 
essential correlations due to the symmetry-breaking mechanism are then 
used as trial states during the projection procedure. The results for ac 
in 80 Zr are consistent with current experimental estimates from the giant 
dipole resonance studies. The preliminary results on the ^-corrections are 
also very encouraging. The calculated values of the nucleus-independent 
Ft = 3069.4(10) and the V ud = 0.97463(24) are consistent with the recent 
evaluations of Ref. [7] . 
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